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Errata 
EM-199 

A  Plane  Wave  Expansion  Theorem  for  Cylindrically 
Radiated  Fields 

a.H.  Karp  and  N.R.  Zitron 


Page  5,  Equation  2(a)  (p  (o)  should  read  (6  (O)) 


Page  7,  line  13 


.(m-p). 


Equation  should  read; 
m-p    X 


XT 


Page  8 


Then  l^f-   should  read : 


ht= 


0 

r  (t  +  1/2) 

2t  +  1 

(kd)  2 


m  odd 

m  even   t  =  0,1,2,, 


Page  9,  Footnote  should  read: 

It  should  he  pointed  out  here  that  the  definition  of  the 
complex  scattering  amplitude  of  the  far  field  f (8 )  in  this 
paper  differs  from  that  of  the  previous  paper  [l].Let 
f  denote  the  complex  scattering  amplitude  of  the  far  field 

in  this  paper  and  f  denote  the  one  in  the  previous  paper. 
The  relation  between  the  two  amplitudes  is  then 


Page  10,  line  7 


Equation  should  read: 


u(r,0)  ~  /— ; —  e 


p    i(kr-  -J-)  ^ 
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Abstract 

An  asymptotic  expansion  for  two-dimensional  outwardly  radiating 
fields  is  developed  from  an  integral  representation  of  these  fields  by 
means  of  a  saddle  point  integration.   The  expansion  is  given  in  terms  of 
inverse  powers  of  the  distance  from  a  point  in  a  fixed  region  A  to  a  point 
in  a  circular  neighborhood  at  a  large  distance  from  A.   The  coefficients 
are  expressed  in  terms  of  plane  waves  and  linear  combinations  of  derivatives 
of  plane  waves  with  respect  to  angle  of  incidence.   The  theorem  may  be 
employed  in  scattering  problems  in  reducing  scattering  of  arbitary  two- 
dimensional  fields  by  arbitrary  cylinders  to  scattering  of  plane  waves  by 
the  arbitrary  cylinders. 


-2- 


Introduction 

The  radiation  from  a  collection  of  cylindrical  sources  within  a 
bounded  cylindrical  domain,  when  viewed  by  an  observer  at  an  infinite  distance, 
cannot  be  distinguished  from  a  plane  wave.   As  the  observer  approaches  closer 
to  the  source  region,  the  nature  of  the  radiation  may  appear  cylindrical  or 
perhaps  even  more  complicated.   If,  however,  the  distance  from  the  source 
is  sufficiently  large,  it  may  be  said  that  the  wave  is  almost,  but  not  quite, 
a  plane  wave. 

This  paper  is  concerned  with  the  relation  between  a  general  cylindrically 
radiated  wave  and  a  plane  wave.   It  is  shown  here  that  the  relation  is 
expressible  as  an  asymptotic  series  involving  plane  waves  and  their  derivatives 
of  various  orders  with  respect  to  angle.   A  few  terms  of  this  series,  obtained 
previously  by  Zitron  and  Karp  [l] ,  were  utilized  to  good  advantage  in  dealing 
with  multiple  scattering  of  plane  waves  by  two  widely  spaced  parallel  cylinders 
of  arbitrary  shape.   The  present  paper  gives  an  expression  for  the  general 
term  of  the  series  for  which  the  coefficient  is  a  linear  combination  of  plane 
waves  and  their  derivatives  with  respect  to  angle  of  incidence.   This 
representation  also  provides  a  more  efficient  way  of  calculating  the  terms 
of  the  expansion  than  was  used  before  [l] .   These  results  were  reported  earlier 
in  abstracted  form  [2]  ,  [3]  . 

The  major  resvilt  presented  in  this  paper  is  an  asymptotic  expansion, 
in  a  distant  finite  region  of  an  integral  representation  for  a  cylindrically 
radiated  field.  Application  of  the  saddle  point  method  to  this  integral  yields 
a  power  series  in  inverse  half- integral  powers  of  kd,  where  k  is  the  propagation 
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Gonstant  and  d  is  a  long  distance  by  which  the  origin  0  of  a  coordinate  system 
located  in  the  source  region  has  been  translated.   Since  the  coefficients  in 
the  asymptotic  expansion  turn  out  to  be  linear  combinations  of  plane  waves 
and  their  derivatives  of  various  orders  with  respect  to  angle  of  incidence, 
the  expansion  can  be  used  to  represent  the  given  field  in  a  new  coordinate 
system  as  a  plane  wave  plus  a  correction  for  the  cujrvature  of  the  wavefront 
in  the  neighborhood  of  a  scatterer  located  near  the  origin  d  of  the  new 
coordinate  system.   The  correction  terms  appear  in  the  form  of  derivatives 
of  plane  waves  with  respect  to  angle  of  incidence.  The  latter  is  merely  a 
parameter.  Thus  the  theorem  has  special  utility  in  scattering  problems,  since 
the  response  to  the  incident  wave  can  be  given  in  terms  of  derivatives  of  the 
response  to  a  plane  wave  as  was  done  by  the  authors  [l]  in  a  previous  paper. 
It  is  possible,  then,  by  means  of  this  expansion,  to  determine  the  scattered 
field  when  radiation  from  a  cylindrical  region  is  incident  upon  a  cylindrical 
scatterer  and  also  to  reduce  multiple  scattering  problems  to  single  scattering 
problems.   Both  of  these  cases  may  be  treated  in  a  quite  general  manner  without 
reference  to  the  detailed  nature  of  the  Incident  radiation  or  the  shape  or 
composition  of  the  scatterer.   In  this  respect,  the  expansion  theorem  may  be 
used  in  place  of  addition  theorems  which  are  applicable  only  when  suited  to 
the  special  geometry  of  the  sources  and  scatterers  and  which,  even  in  such  cases, 
may  be  clumsier  to  apply. 
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Derivation  of  the  Expansion 


The  geometry  of  the  situation  is  illustrated  in  Fif^ire  1. 


Figure  1 

Let  A  be  a  fixed  region  of  arbitrary  size  and  let  d  be  the  distance 
from  a  point  0  at  the  center  of  the  smallest  circle  containing  A  to  the  center 
0'  of  a  circle  B.   Let  the  diameters  of  the  circles  A  and  B  be  small  in  comparison 
with  d.   Let  P  be  a  point  in  B  other  than  0'.   Let  r  andg  be  the  polar 
coordinates  of  p  with  origin  at  0^  where  9   is  measured  from  the  0-0*  axis. 
Let  X  and  y  be  cartesian  coordinates  of  p  with  origin  at  0'.   The  principal 
result  of  this  paper  may  then  be  stated  as  follows: 

Theorem  Consider  any  two-dimensional  radiated  field*  represented  by 


*  A  radiated  field  is  a  solution  of  the  reduced  wave  equation  which  satisfies 
the  Sommerfeld  radiation  condition  for  an  outgoing  wave. 
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(1) 


u(r,0) 


f(p)e^^  cos(0-p)  ^p 


in  the  exterior  of  a  bounded  domain  A  with  center  0  containing  any  combination 
of  sources  and  scatterers  where  C  is  the  Sommerfeld  contour  for  h!;   (kr)  and 
f(p)  is  analytic.   Then  u(r,0)  has  an  asymptotic  expansion  in  any  neighborhood 
P  of  a  point  0'  such  that  p  is  disjoint  with  respect  to  the  smallest  circle 
containing  A  provided  that  kd  >  >  1  where  d  is  the  distance  between  0  and  0' 
(Figtire  1).   The  expansion  is  given  in  terms  of  plane  waves  and  their  derivatives 
with  respect  to  angle  as  follows: 


(2) 


u  ^  e 


ikd 


t=o  "^^-    (kd) 


t+±/2       ^°^  ^"y  ^' 


where 


,(2t), 


2t   2t-p   X 


(2a)  h^^^''(o)=y-  V-"  V- 

p=0  X=0   j=0 


(^^-P)(B  (0))(P)  f(^0)'^  M 


and 


r(t)  is  the  Gamm  function,  c'^  =   ,V   .,  ,   v(p)  =  e^^^^  =°"P"^  "^^^^ 


f(0) 


,^-J 


Lip' 


^-j 


f(P) 


p=o 


,  (^3(0)) 


(p) 


,p+l 


ds 


P+1 


U(s) 


s=0. 


=(P^o))(2^-P^ 


2t-p 


ds 


2t-p 


P"(s) 


s=0 


P(s)  =arcos  (l+  is^)  ,  &(0)   =  ff 


ds 


s=0 


Proof 

Let 


u(r,0)  =J      e'^=°^^^-P^  f(p)dp 


^1 
where  f (p)  is  analytic  and  C^  is  the  Soramerfeld  contour  employed  in  the  integral 
representation  [5] , [6]  of  H^  \kr)  (See  figure  2),   JThe  time  dependence  is 
assumed  to  be  e"''"  J    Expansion  of  cos(0-P)  in  the  exponent  of  (l)  in  terms 
of  X,   J,   and  d  (figure  l)  yields 

(3)     u(r,e)  =f      e^^^^  '^^^P-^  ^^"P^  f(p)e^^^  ^^^^dp  =  f    s{^)e'^'''°'% 
^1  ^'l 

where  gO)  =  e^^^^  "°^P-*^  "'^^^fO)  * 

It  is  clear  that  e'^  c^sp  .j^^^^^   rapidly  for  large d  and  that  the 
saddle  point  method  is  appropriate  for  evaluation  of  the  integral.   No 

difficulty  is  encountered  in  shifting  the  contoior  as  d  becomes  large  since 

2 
the  integrand  is  analytic.   The  transformation  cos  p  =  1+  is  transforms  C^ 

into  a  steepest  descent  contour  (figure  2)  and  (3)  into 

(k)         u  =  e^^<^f  g(p(s))  If    e-^'^^'  ds  =  e^^^  [   h^s)]  e  "^^^^  ds 


Sommerfeld  Contour       Saddle  Contour  of  Steepest  Descent 


Figure 


-7- 


vhere  h[p(s)]  .  e^^^^  ^^^^^^^"^  ''^   ^^^^^fO(s))  ff]    ' 

If  h  f3(s)  is  expanded  in  a  Taylor  series  about  t"      lie  point 
s  =  0  and  a  finite  number  of  these  terms  are  retained,  it  is  clear  that 


ikd 


(5)  u^e^^^  Yi 


m=0 


h 


(in 


k 


m: 


m  -kds 
s  e      ds 


IM^  |/'")(0) 


m=0 


m 


m 


where  I  =  /   s"^  e"    ds  and  the  superscribed  parenthesis  e" 

loo 

differentiation  with  respect  to  s. 

The  function  h   (s)  can  be  represented  in  terms  of  1 
plane  waves  as  follows: 


:ives  of 


h 


(m)  ^  ^^^  ^(m)  ^  Y^   ^m  ^(m-p)   (p) 

p  '^       s 


P=0 


by  Leibnitz's  Rule  [t]  . 

A  rule  given  by  Qoursat  [8]  for  higher  order  derivatives  of  implicit 

functions  yields 


X=0   e=0 


X! 


where  the  square  bracketed  superscript  signi 


fies  differentiation  with  respect 


to  p.   But  p(0)  =  0.   Therefore, 


.        .               m-£        (P   (0))'      ^  r^^i 

(ni-P)(0)=f2       Xi ^        ^°^ 

x=o 


A  further  application  of  Leibnitz's  rule  yields 


3=0 


Therefore 


h^^^(O) 


m       m-p  X 

p=0     X=0        j=0 


c"^c^ 


-^  (p'(o))^-p^  p^p^(o)  f\ly'^  vW  • 


Evaluation  of  I     results   in 

m 


m  odd 


(     r(2li) 

m  even 

m+l 

V.(kd)  ^ 

Let  m  =  2t 

r° 

m  odd 

Then  !_,       =   ■" 
2t 

1 

r(t->|) 

m+l 
2 

m  even 

(M) 


t   =  0,1,2, 


and 


(6) 


u^e 


iM     "       h(gt)(o)       r(t4) 


Z  T 


2tTy 


t=0 


(kd) 


t+1/^ 


where 


(7 


2t        a:-p       X         C^  C^ 


)      h^^)(o)=  f;^  Y^  fl   ^ir^  (p'^°))^^"''^(Ps^°))^^^f^'''^(o)vtj^co) 


p=0     X=0        j=0 


It   is  convenient  to  note  that 
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(8)  r(p) 


P^^'(O) 


0     p  odd 


*>  2  -^2 


p: 


.P/P 


2^(^)V 


i  (2v)j\    ^ 

e 

2v   , 

p  even 

where  p  =  2V , 

V  =  0,1,2 

The  symbol  p  (O)     ,   evaluated  in  (8)  means  the  p   derivative  of   ■^■^  ' 


evaluated  at  s  =0. 

It  is  clear,  then,  from  (2)  and  (2a)  that 


?ients  in  the 


form  of  derivatives  of  plane  waves  with  respect  to  angle  obtained  previously  [l] 
by  the  authors  were  not  accidental. 


It  should  be  pointed  out  here  that  the  definition  of  the  complex  scattering 
amplitude  of  the  far  field  f (p)  in  this  paper  differs  from  that  of  the  previous 
paper. L^J   Let  f  denote  the  complex  scattering  amplitude  of  the  far  field  in 
this  paper  and  f  denote  the  one  in  the  previous  paper.   The  relation  between 
the  two  amplitudes  is  then  .jt 
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Application  of  the  Theorem,  and  Discussion. 

We  shall  now  give  an  illustrative  calculation  in  outline  form. 
Comments  on  the  nature  of  the  work  are  also  included.   For  an  extended 
illustration  see  [4]  . 

1.   The  function  f(p)  in  the  integral  representation  (l)  for  u  is  related 
to  the  scattering  amplitude  by  the  formulas. 

i-r-   ikr—  ^ 


rtkr 


i  f(9)  Hf(0) 


2. 


d^ 
ds 


(s) 


72   e"^ 


s=0 


3.   To  apply  the  formula  (2a)  to  a  radiated  field  it  is  sufficient  to 
know  f(0)  which  determines  the  function  u.   Thus  the  radiated  field  need 
not  be  given  in  the  form  of  our  contour  integral  in  order  to  apply  our  result. 
Instead,  it  suffices  to  know  f (0)  in  a  small  neighborhood  of  0  =  0. 

It  is  not  necessary  to  know  f (s)  in  a  large  angular  region.   In  this 
respect  again  our  result  differs  from  the  more  usual  expansion  theorem,  which 
proceeds  in  inverse  powers  of  r. 

h.      It  is  convenient  to  bear  in  mind  that  the  functions  v"-  ^  (O)  really 
contain  the  new  coordinates  x  and  y.   For  example 

,^,    ik(x  cos6  +  y  sinS) 
v(Oj  =  e        .    J    I-/ 
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5. 


v(0)  =  e 
vW  (0) 


ikx 


ikye 


ikx  _  dv(p)| 


^^  P=0  / 


vl2J(0)  =  ik[-x  +  iky^J 


2t  ikx   d  V 
e    =  — 

dp 


3=0. 


,  etc. 


A), 


Consider  for  example  the  expanGion  of  the  function  H   (kr). 


We  have 


i(kr4)" 


nf  ^(^^  -  J^     ^ 


1+  0(i) 
r 


r— >  00 


Hence  f(0)  h  1,      f (e)  =  -  .   Applying  (2a)  we  find 


l(M-f) 


ikx  e        +  Jiik[-x-Hky]e  ^. 


2_   i(^<^-!> 

jtkd     ^ 


v(0)      +     |rj^^    (y(0)    +  W'-2J(o)) 


+  Uv^^^fO))    + 


We  shall  now  explain,  in  brief,  how  the  expansion  theorem  may  be  employed 
in  calciilating  the  diffraction  of  an  arbitrary  radiating  field. 

In  such  a  calculation  the  function  f (S)  for  the  incident  radiation  will 
be  known.  The  incident  radiation  will  then  be  represented  in  the  form  (2a), 
in  a  coordinate  system  x,y  fixed  in  the  neighborhood  of  the  scatterer. 

In  addition  we  must  know  the  response  amplitude  f  {d,9   )  say,  of  the 
scatterer,  in  direction  9,   due  to  excitation  by  an  incident  plane  wave 
exp  ik[_x  COS©  +  y  sin0  1  =  v(e  ).   Since  0  is  a  parameter,  we  can  then 
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calculate  the  response  to  any  of  the  functions  v*-  -'  (O)  in  termsof  which  the 
incident  radiation  will  have  been  expressed,  since  the  response  to  a  derivative 
of  a  plane  wave  is  the  derivative  of  the  response  to  the  plane  wave,  each 
term  of  the  response  to  this  incident  radiation  will  have  a  coefficient 
identical  to  the  one  which  appears  in  the  expansion  theorem  (a  linear  combination 
of  derivatives  with  respect  to  angle  of  incidence)  operating  upon  the  response 
to  a  plane  wave  of  a  given  angle  of  incidence.   For  a  detailed  illustration 
cf[ij-J  and  [l]. 
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